Transitional Lu and Spherical Ta Ground-State Proton Emitters 

in the Relativistic Hartree-Bogoliubov model 

G.A. Lalazissis^, D. Vretenar^'^, and P. Ring^ 
^ Physik- Department der Technischen Universitdt Miinchen, D-8574-8 Garching, Germany 
^ Physics Department, Faculty of Science, University of Zagreb, Croatia 



"q>^' ' Properties of transitional Lu and spherical Ta ground-state proton emit- 

On ' ters are calculated with the Relativistic Hartree Bogoliubov (RHB) model. 
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The NL3 effective interaction is used in the mean-field Lagrangian, and pair- 
ing correlations are described by the pairing part of the finite range Gogny 
interaction DIS. Proton separation energies, ground-state quadrupole defor- 
mations, single-particle orbitals occupied by the odd valence proton, and the 
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■ corresponding spectroscopic factors are compared with recent experimental 

o ■ 

On ■ data, and with results of the macroscopic-microscopic mass model. 
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The structure and decays modes of nuclei beyond the proton drip-line represent one of the 
most active areas of experimental and theoretical studies of exotic nuclei with extreme isospin 
values. In the last few years many new data on ground-state and isomeric proton radioactiv- 
ity have been reported. In particular, detailed studies of odd-Z ground-state proton emitters 
in the regions 51<Z<55 and 69<Z< 83, have shown that the systematics of spectroscopic 
factors is consistent with half-lives calculated in the spherical WKB or distorted-wave Born 
(DWBA) approximations More recent data indicate that the missing region of 

light rare-earth nuclei contains strongly deformed systems at the drip-lines. 

In the theoretical description of ground-state and isomeric proton radioactivity, two es- 
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sentially complementary approaches have been reported. One possibihty is to start from a 
spherical or deformed phenomenological single-particle potential, a Woods-Saxon potential 
for instance, and to adjust the parameters of the potential well in order to reproduce the ex- 
perimental one-proton separation energy. The width of the single-particle resonance is then 
determined by the probability of tunneling through the Coulomb and centrifugal barriers. 
Since the probability strongly depends on the valence proton energy and on its angular mo- 
mentum, the calculated half-lives provide direct information about the spherical or deformed 
orbital occupied by the odd proton. For a spherical proton emitter it is relatively simple 
to calculate half-lives in the WKB or DWBA approximations [0]. On the other hand, it 
is much more difficult to quantitatively describe the process of three-dimensional quantum 
mechanical tunneling for deformed proton emitters. Modern reliable models for calculat- 
ing proton emission rates from deformed nuclei have been developed only recently A 
shortcoming of this approach is that it does not predict proton separation energies, i.e. the 
models do not predict which nuclei are likely to be proton emitters. In fact, if they are used 
to calculate decay rates for proton emission from excited states, the depth of the central 
potential has to to be adjusted for each proton orbital separately. In addition, the models 
of Refs. 0,^,01 do not provide any information about the spectroscopic factors of the proton 
orbitals. Instead, experimental spectroscopic factors are defined as ratios of calculated and 
measured half-lives, and the deviation from unity is attributed to nuclear structure effects. 

In Refs. we have used the relativistic Hartree Bogoliubov (RHB) theory to calculate 
properties of proton-rich spherical even-even nuclei with 14<Z<28, and to describe odd-Z 
deformed ground-state proton emitters in the region 53 < Z < 69. RHB presents a relativis- 
tic extension of the Hartree-Fock-Bogoliubov theory, and it provides a unified framework 
for the description of relativistic mean- field and pairing correlations. Such a unified and self- 
consistent formulation is especially important in applications to drip-line nuclei. The RHB 
framework has been used to study the location of the proton drip-line, the ground-state 
quadrupole deformations and one-proton separation energies at and beyond the drip line, 
the deformed single particle orbitals occupied by the odd valence proton, and the corre- 
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spending spectroscopic factors. The results of fully self-consistent calculations have been 
compared with experimental data on ground-state proton emitters. However, since it is very 
difficult to use the self-consistent ground-state wave functions in the calculation of proton 
emission rates, one could say that the RHB model provides informations which are com- 
plementary to those obtained with the models of Refs. It should be noted that in 
the relativistic framework the strength and the shape of the spin-orbit term are determined 
selfconsistently. This is essential for a correct description of spin-orbit splittings in regions 
of nuclei far from stability, where the extrapolation of effective strength parameters becomes 
questionable. The motivation for the present work are the very recent data on proton emis- 
sion from the closed neutron shell nucleus ^^^Ta 0, and the proposed experiment to search 
for direct proton emission from ^^^Lu |]10[- The analysis of ground-state proton radioactivity 
in the Lu and Ta isotopes completes our study of deformed and transitional proton emitters 
in the region 53 < Z < 73. 

A very detailed description of the relativistic Hartree-Bogoliubov theory can be found, 
for instance, in Ref . ^ . In the following we only outline the essential features of the model 
that will be used to describe nuclei at the proton drip-line. The ground state of a nucleus is 
represented by the Slater determinant of independent single-quasiparticle states, which are 
obtained as solutions of the relativistic Hartree-Bogoliubov equations 
fiD-m-X A \ /f/fc(r)\ /f/fc(r)\ 

-A* -ho + m + \) [Vkir)) \Vkir))' 
The column vectors denote the quasi-particle spinors and Ek are the quasi-particle energies. 
In the Hartree approximation for the self-consistent mean field, the single-nucleon Dirac 
Hamiltonian reads 

h^ = -ta-V + (3{m + g^a{r)) + g^T^u\v) + g^p^v) + e ^^ ~ ^'^ A°(r). (2) 

It describes the motion of independent Dirac nucleons in the mean-field potentials: the 
isoscalar scalar cr-meson potential, the isoscalar vector cj-meson, and the isovector vec- 
tor p-meson potential. The photon field A accounts for the electromagnetic interaction. 
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The meson potentials are determined self-consistently by the solutions of the corresponding 
Klein-Gordon equations. The source terms for these equations are calculated in the no-sea 
approximation. Because of charge conservation only the third component of the isovector 
p-meson contributes. For an even-even system, due to time reversal invariance the spatial 
vector components a;, ps and A of the vector meson fields vanish. In nuclei with odd num- 
bers of protons or neutrons time reversal symmetry is broken. The odd particle induces 
polarization currents and the time-odd components in the meson fields. These components 
play an essential role in the description of magnetic moments and of moments of inertia in 
rotating nuclei. However, their effect on deformations and binding energies is very small and 
can be neglected to a good approximation. As in our previous studies of nuclei at the proton 



drip-lines, we choose the NL3 set of meson masses and meson-nucleon coupling constants |11 



for the effective interaction in the particle-hole channel: m = 939 MeV, mg. = 508.194 MeV, 
= 782.501 MeV, rrip = 763.0 MeV, = 10.217, g2 = -10.431 fm-^, ^3 = -28.885, 
= 12.868 and gp = 4.474. 
The pairing field A is defined 

Aab{r, r') = ^Y. ^abcdir, r')Kcd(r, r'), (3) 

c,d 

where Vabcdii^^, r') are matrix elements of a two-body pairing interaction, and Kcd(r, r') is the 
pairing tensor. The pairing part of the phenomenological Gogny force 

VPP{1,2) = e-«"^~"^)/^»)' {Wi + BiP" - HiP^ - MiP^P^), (4) 

1=1.2 



with the set DIS [T^ for the parameters /Xj, Wi, Bi, Hi and Mi {i = 1, 2), is used to describe 
pairing correlations. 

The RHB single-quasiparticle equations (|1]) are solved self-consistently. The iteration 
procedure is performed in the quasi-particle basis. The chemical potential A has to be 
determined by the particle number subsidiary condition in order that the expectation value 
of the particle number operator in the ground state equals the number of nucleons. A simple 
blocking prescription is used in the calculation of odd-proton and/or odd-neutron systems. 



The blocking calculations are performed without breaking the time-reversal symmetry. The 
resulting eigenspectrum is transformed into the canonical basis of single-particle states, 
in which the RHB ground-state takes the BCS form. The transformation determines the 
energies and occupation probabilities of the canonical states. 
The one-proton separation energies 

Sp{Z,N) = B{Z,N) - B{Z -1,N) (5) 

for the Lu and Ta isotopes are displayed in Fig. |l], as function of the number of neutrons. 
The predicted drip-line nuclei are ^^'^Lu and ^^^Ta. In the process of proton emission the 
valence particle tunnels through the Coulomb and centrifugal barriers, and the decay prob- 
ability depends strongly on the energy of the proton and on its angular momentum. In 
rare-earth nuclei the decay of the ground-state by direct proton emission competes with f3~^ 
decay; for heavy nuclei also fission or a decay can be favored. In general, ground-state pro- 
ton emission is not observed immediately after the drip-line. For small values of the proton 
separation energies, the width is dominated by the (3~^ decay. On the other hand, large 
separation energies result in extremely short proton-emission half-lives, which are difficult 
to observe experimentally. For a typical rare-earth nucleus the separation energy window 
in which ground-state proton decay can be directly observed is about 0.8 - 1.7 MeV p|. In 
Fig. ID we have compared the calculated separation energies with experimental transition 
energies for ground-state proton emission in ^^''Lu, ^^^Lu |T^, ^^^Ta 0, ^^^Ta |]14|, and ^^''Ta 



p!5| . In all five cases an excellent agreement is observed between model predictions and ex- 
perimental data. In addition to ^^^Lu, which was the first ground-state proton emitter to be 



discovered |T6|, and ^^''Lu, the self-consistent RHB calculation predicts ground-state proton 



decay in ^"^^Lu. The calculated one-proton separation energy —1.77 MeV corresponds to a 
half-life of a few fis, if one assumes that the nucleus is spherical. Direct proton emission 
with a half-life of the order of few /is is just above the lower limit of observation of cur- 
rent experimental facilities. An experiment to search for direct proton emission from ^^^Lu 
has been proposed recently [|lOl- For the Lu ground-state proton emitters, in Table | the 



results of the RHB model calculation are compared with the predictions of the finite-range 
droplet (FRDM) mass model: the projection of the odd-proton angular momentum on the 
symmetry axis and the parity of the odd-proton state Qp ||T^, the one-proton separation 
energy |]T^, and the ground-state quadrupole deformation |[T^. We have also included the 



RHB spectroscopic factors, and compared the separation energies with the experimental 
transition energies in ^^°Lu and ^^^Lu. Both theoretical models predict oblate shapes for 
the Lu proton emitters, and similar values for the ground-state quadrupole deformations. 
On the other hand, while the FRDM assigns spin and parity 5/2^ to the deformed single- 
particle orbitals occupied by the odd valence proton in all three proton emitters, the RHB 
model predicts the 7/2^ [523] Nilsson orbital to be occupied by the odd proton. We also 
notice that the RHB separation energies are much closer to the experimental values. The 
spectroscopic factors of the 7/2^ [523] orbital are displayed in the sixth column of Table |. 
The spectroscopic factor of the deformed odd-proton orbital k is defined as the probability 
that this state is found empty in the daughter nucleus with even number of protons. 

In the detailed analysis of odd-Z proton emitters (53 < Z < 69) [0] it has been shown 
that, while the proton-rich isotopes of La, Pr, Pm, Eu and Tb are all strongly prolate de- 
formed {P2 ~ 0.30 — 0.35), Ho and Tm isotopes at the proton drip-line display a transition 
from prolate to oblate shapes. Spherical shapes are expected as the nuclei with unbound 
protons approach the = 82 neutron shell. The Lu proton emitters are found in the tran- 
sitional region between oblate and spherical shapes. This is illustrated in Fig. where we 
plot the binding energy curve for ^^^Lu as function of the quadrupole deformation param- 
eter. The binding energies result from self-consistent RHB/NL3 calculations performed by 
imposing a quadratic constraint on the quadrupole moment. A very shallow minimum is 
found at /3 ~ —0.15, but otherwise the potential is rather flat with a shoulder at /5 = 0. 
In Fig. ^ we compare the ground-state quadrupole deformations of the proton-rich Lu iso- 
topes with those of Ho and Tm ^ . For < 80 all three chains of isotopes display oblate 
deformations; starting with A^ = 81 a sharp transition to the spherical shape is observed. 

The proton-rich Ta isotopes are spherical. In Fig. |l| we compare the calculated one- 
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proton separation energies with experimental transition energies for Ta M, Ta [14 



and ^^''Ta |T^. The predictions for the spherical orbitals occupied by the odd proton and 
the corresponding spectroscopic factors are displayed in Table ^ Results of the FRDM cal- 
culation Jl^ have also been included in the comparison. As in the case of the Lu ground-state 
proton emitters, an excellent agreement between RHB separation energies and experimental 
data on transition energies for proton emission is observed. In particular, our calculation 
reproduces the very recent data on proton emission from the closed neutron shell nucleus 
^^^Ta 1^. The significant decrease in proton binding for ^^^Ta, as compared to ^^^'^^^Ta, has 
been associated with the = 82 closure. In comparison, the FRDM results are found to be 
in rather poor agreement with experimental data. Except for ^^^Ta, the spherical orbitals 
predicted to be occupied by the odd proton agree with the experimental assignments, and 
the theoretical spectroscopic factor of the /iii/2 orbital in ^^^Ta is very close to the exper- 
imental value. For ^^^Ta the RHB model predicts ground-state proton emission from the 
^11/2 orbital. The experimental assignment for the ground-state configuration is Si/2, but 
an alpha decaying state is identified in ^^'^Ta at an excitation energy of only 22(5) keV and 
assigned to an /iii/2 isomer |TB[. We have also calculated the one-proton separation energy 
for ^^^Ta™": Sp = —1.250 MeV, the orbital is /iii/2 and the spectroscopic factor is 0.79. 



This is to be compared with the experimental transition energy Ep = 1.103(12) MeV [119 
assigned to the /iii/2 orbital with the experimental spectroscopic factor 0.92(4) [jl||. 

In conclusion, the relativistic Hartree-Bogoliubov model has been applied in the descrip- 
tion of ground-state properties of transitional Lu and spherical Ta proton emitters. The 
NL3 effective interaction has been used for the mean-field Lagrangian, and pairing correla- 
tions have been described by the pairing part of the finite range Gogny interaction DIS. We 
would like to emphasize that this particular combination of effective forces has been used 
in most of our recent applications of the RHB model, not only for spherical and deformed 
/9-stable nuclei, but also for nuclear systems with large isospin values on both sides of the 
valley of /5-stability. The model parameters therefore have not been adjusted to the specific 
properties of nuclei studied in this work, or to the properties of deformed proton emitters 



discussed in Refs. The self-consistent calculation reproduces in detail the observed 

transition energies for ground-state proton emission in ^^°Lu, ^^^Lu |jl3[, ^^^Ta 0, ^^^Ta 
I^, and ^^^Ta as well as the assignments for the orbitals occupied by the valence odd 



proton and the corresponding spectroscopic factors. The model also predicts the one-proton 
separation energy of —1.77 MeV for the possible proton emitter ^^^Lu. Oblate ground-state 
deformations are predicted for all Lu proton emitters, while spherical shapes are calculated 
for the Ta isotopes at and beyond the proton drip-line. With the excellent agreement ob- 
served between the RHB results and the very recent data on proton emission from the closed 
neutron shell nucleus ^^^Ta P], we complete our study of ground-state properties of deformed 
and transitional proton emitters [^|[. 
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Figure Captions 

FIG. 1. Proton separation energies for Lu and Ta isotopes at and beyond the drip-line. Re- 
sults of self-consistent RHB calculations are compared with experimental transition energies for 



ground-state proton emission in ^^"Lu, ^^^Lu |13|, ^^^Ta M, ^^^Ta |14|, and ^^^Ta |15]. Filled 



diamonds denote the negative values of the transition energies Ep. 

FIG. 2. Potential energy curve for ^^^Lu. The self-consistent RHB calculations are performed 
with constrained quadrupole deformation. 

FIG. 3. Self-consistent ground-state quadrupole deformations for Ho, Tm and Lu isotopes, at 
and beyond the proton drip-line. 
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TABLE I. Lu ground-state proton emitters. Results of the RHB calculation for the one-proton 
separation energies Sp, quadrupole deformations /?2, and the deformed single-particle orbitals oc- 
cupied by the odd valence proton, are compared with predictions of the macroscopic-microscopic 
mass model, and with the experimental transition energies. All energies are in units of MeV; the 
RHB spectroscopic factors are displayed in the sixth column. 



N 



Sn 



P2 



p-orbital 



Sp P2 Ep exp. 



149Lu 78 
isoLu 79 
is^Lu 80 



-1.77 -0.158 7/2- [523] 0.60 5/2" -1.51 -0.175 

-1.31 -0.153 7/2- [523] 0.61 5/2" -1.00 -0.158 1.261(4) |l| 

-1.24 -0.151 7/2- [523] 0.58 7/2" -0.99 -0.150 1.233(3) H 



TABLE II. Spherical Ta ground-state proton emitters. RHB results for the proton separation 
energies, the single-particle orbitals occupied by the odd proton, and the corresponding spectro- 
scopic factors are compared with the predictions of the finite-range droplet (FRDM) mass model 
and with experimental data. 



RHB/NL3 


i55Ta 




-1.677 






h 11/2- 




spectr. factor 


0.60 




Sp 


-1.129 






d 3/2+ 




spectr. factor 


0.51 


i57Ta 


Sp 


-1.040 






h 11/2- 




spectr. factor 


0.42 



FRDM |T7| 



EXP. 



-1.09 

9/2- 

-0.60 

3/2" 

-0.48 

9/2- 



-1.765(10) i 
11/2- 
0.58(20) 
-1.007(5) 

3/2+ 

-0.927(7) m 

1/2+ 
0.56(24) 
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